function output = f(lambda)

p0 = -6;

output = lambda - 9.8*tanh(abs(p0)/(lambda^.5));

end
